Phonon Dispersion Relations in PrBa 2 Cu 3 06+x (x ~ 0.2) 
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We report measurements of the phonon dispersion relations in non-superconducting, oxygen- 
deficient PrBa2Cu306+a: (x ~ 0.2) by inelastic neutron scattering. The data are compared with a 
model of the lattice dynamics based on a common interatomic potential. Good agreement is achieved 
for all but two phonon branches, which are significantly softer than predicted. These modes are 
found to arise predominantly from motion of the oxygen ions in the CUO2 planes. Analogous modes 
in YBa2Cu306 are well described by the common interatomic potential model. 

PACS numbers: 74.72.-h, 63.20.Dj, 61.12.Ex 



I. INTRODUCTION 
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The suppression of superconductivity by 
praseodymium in PrBa2Cu306+x has for some years 
presented an outstanding problem in the field of cuprate 
superconductivity^^ Although superconductivity at 
T c w 90 K is exhibited by the majority of compounds 
with the composition i?Ba2Cu30y, where R is a 
rare-earth, PrBa2Cu306+x combines antiferromagnetic 
ordering with semiconducting resistivity across the 
whole of the known phase diagram (0 < x < 1). 

An influential model of the electronic structure of 
PrBa 2 Cu306+^ suggests that the absence of supercon- 
ductivity is caused by localisation of holes due to hybridi- 
sation of the Pr 4/ and O 2p orbitals. This is supported 
by the observation of an enhanced ordering temperature 
Tp r of the Pr sublattice, which is an order of magnitude 
larger than the ordering temperature Tr of the rare-earth 
sublattice in other i?Ba2Cu307 compounds. 

Neutron scattering studies of the magnetic structure* 
and dispersive magnetic excitations 5,6 have probed 
the strengths of the different magnetic couplings in 
PrBa 2 Cu306+ a: , and indicate that the high value of Tp r is 
partly due to an enhanced Pr-O-Pr superexchange and 
partly due to the presence of a significant coupling be- 
tween the Cu and Pr sublattices. Both of these effects 
could be the result of hybridisation between Pr 4/ elec- 
trons and the Cu02 planes. 

Hybridisation is expected to affect the phonon disper- 
sion relations through the Pr-0 bonds. In this paper we 
present measurements of the phonon dispersion relations, 
performed by inelastic neutron scattering on an oxygen- 
deficient single crystal of PrBa2Cu306+ K (x « 0.2). The 
data are compared with a model of the lattice dynam- 
ics based on a common interatomic potential, which was 
originally developed for YBa2Cu30g and later modified 



II. EXPERIMENTAL DETAILS 

The single crystal sample of PrBa2Cu306+z was pre- 
pared by top seeding a flux and had a mass of ~ 2 g. 
It was reduced at 700° C in a flow of 99.998% argon for 
100 hours and quenched to room temperature resulting 
in an estimated oxygen content of x rs 0.2. Using x-ray 
Laue diffraction it was aligned with the [110] direction 
vertical. It was then glued onto an aluminium mount 
using G.E. varnish and held securely in place with alu- 
minium wire. Preliminary neutron diffraction measure- 
ments showed the crystal mosaic to be ~ 1°. 

The neutron scattering experiments were performed on 
the INI and IN22 triple-axis spectrometers at the Insti- 
tut Laue-Langevin research reactor. The configuration of 
INI was as follows: flat copper (200) monochromator, py- 
rolytic graphite (002) analyser with horizontal and slight 
vertical focussing, no collimation before the monochro- 
mator, 60' collimator between the monochromator and 
sample, and pyrolytic graphite filter between the sam- 
ple and analyser to reduce contamination of the beam by 
second and higher order reflections from the monochro- 
mator. The spectrometer was operated in the constant 
fcf configuration with a final energy of 34.8 meV. Low 
efficiency beam monitors, placed before the monochro- 
mator and between the filter and the analyser, were used 
to gauge the incident beam intensity and to check for 
accidental Bragg scattering. The crystal was mounted 
inside a displex refrigerator and aligned such that the 
[110] and [001] directions lay within the scattering plane. 
All measurements were performed at T = 12 K. The con- 
figuration of IN22 was identical except for the use of a 
vertically curved pyrolytic graphite (002) monochroma- 
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FIG. 1: Plots (a), (b), (c) and (d) show typical energy scans. The points are the raw neutron scattering data (averaged over 
several runs with different, overlapping energy ranges), and the lines are Gaussian fits. The value of Q is shown at the top 
of each plot. In plots (a) and (b), from INI, a monitor count of 2 x 10 6 corresponds to a counting time of 2-5 minutes per 
point (depending on energy), while in plots (c) and (d), from IN22, a monitor count of 500 corresponds to ~ 2 minutes per 
point. Plots (e), (f), (g) and (h) show simulations of the energy-scans, produced from the common interatomic potential model 
described in the text, which facilitate the assignment of symmetry labels to the different phonon branches. 



tor, no collimators and a fixed final energy of 14.7 me V. A 
low efficiency beam monitor placed before the monochro- 
mator was used to gauge the incident beam intensity, but 
there was no monitor between the filter and the analyser 
to check for accidental Bragg scattering. The crystal was 
mounted inside a helium cryostat and aligned such that 
the [110] and [001] directions lay within the scattering 
plane. Most of the measurements were performed at a 
temperature of 1.6 K. 



III. MEASUREMENTS 

The phonon dispersion relations in PrBa2Cu306.2 were 
measured by performing a series of constant Q en- 
ergy scans at positions across the Brillouin zone (Q 
is the neutron scattering vector). On INI, scans 
were performed over an energy range fuo = 30-90 meV 
with Q = (3±/i, 3±/i, 0) to excite predominantly lon- 
gitudinal modes and Q = (h, h, 15) to excite predom- 
inantly transverse modes. On IN22, scans were per- 
formed over an energy range Tuo = 0-24 me V with 
Q = (2 — h, 2—h, 0), (0, 0, 7+1) to excite predominantly 
longitudinal modes and Q = (h, h, 7), (h, h, 8), (2, 2, I) 
to excite predominantly transverse modes. Large values 
of Q were chosen to maximise the phonon scattering cross 
section (which is proportional to Q 2 ), while minimising 
the cross section for magnetic crystal field excitations 
(which scales with the magnetic form factor, falling off 
at high Q)& 



The scattering intensity of a phonon mode is deter- 
mined by its dynamical structure factor, a quantity that 
varies within each Brillouin zone and also between differ- 
ent zones. This means that some modes have a measure- 
able intensity only at certain points within each Brillouin 
zone. In such cases, or when a scan was contaminated 
by accidental Bragg scattering, energy scans were per- 
formed in neighbouring (±h) Brillouin zones to gain as 
many points on the dispersion curves as possible. 

Some typical scans taken on INI and IN22 are shown 
in Figs.[I](a)-(d). The peaks are quite broad because of 
the horizontal focussing of the analysers. There is little 
evidence of magnetic scattering due to crystal field exci- 
tations, which are known to exist in PrBa2Cu306 at en- 
ergies of 61.5, 65.2, 76.0 and 84.7meV;S, but it is possible 
that a small amount of magnetic scattering might con- 
tribute to the broadness of the observed phonon branches 
in the energy range 70-90 meV. 

To determine the dispersion of the observed phonon 
modes in an unbiased way, each energy scan was fitted 
with a lineshape constructed from several gaussians su- 
perimposed on a constant background. The centres, am- 
plitudes and linewidths of the gaussians and the height 
of the background were refined using a least squares 
method. Figure [21 shows the dispersion curves obtained 
by plotting the peak centres as a function of the phonon 
wavevector q, where q is obtained from 



Q = q + T, (l) 
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and r is a reciprocal lattice vector. 
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FIG. 2: Phonon dispersion curves in PrBa2Cu306.2- The 
points are the fitted peak centres obtained from the constant 
Q energy scans. The lines are the dispersion curves predicted 
by the common interatomic potential model, adapted from 
YBa2Cu3C>6 to PrBa2Cu3C>6. The x-axis is labelled with 
the h or I value of q, where q = (hhO) or (00Z) , in recip- 
rocal lattice units (rlu). The Ei modes were measured at 
Q = (3±h, 3±ft, 0), {h, h, 15) and (2 -ft, 2-h, 0). The E 4 
modes were measured at Q = (h, h, 7) and (h, h, 8). The Ai 
and A 5 modes were measured at Q = (0, 0, 7+0 and (2, 2, I) 
respectively. Branches Ei(l) and Ei(4) are plotted as open 
circles and squares respectively. Both are significantly shifted 
with respect to the predictions of the model. The open trian- 
gles in the Ai plot indicate modes that could not be indexed 
by comparison with the model. On IN22 there was no moni- 
tor to check for accidental Bragg scattering, so it is possible 
that this was the cause of these peaks. 



IV. MODEL 

A model based on a common interatomic potential was 
previously developed for the phonon dispersion curves in 
YBa 2 Cu306>£ We adapted this for PrBa2Cu306 by in- 
serting the appropriate lattice parameters, atomic masses 
and nuclear scattering lengths, and the model was used 
to calculate the frequencies and dynamical structure fac- 
tors of the phonon modes at the points in reciprocal space 
where our measurements were performed. 

Comparison of the Q-variation of the predicted and 
observed modes allowed us to identify the symmetries of 
the different branches. In this process the predicted fre- 
quencies were used only as a guide. More weight was 
given to the correspondence between the predicted and 
observed intensities of the modes. To facilitate the com- 
parison, simulations of the energy-scans were produced 
by combining the predicted dynamical structure factors 
with the spectrometer resolution function, calculated us- 
ing the approximate Cooper-Nathans method^ 

Figs, n ( e ) - (h) show the simulations corresponding to 
the selected raw data shown in plots (a)-(d) . Most of the 
observed modes were in good qualitative agreement with 
the simulations, so assignment of the symmetries was 
clear. However, plots (b) and (f) show a clear discrepancy 
between the observed and predicted phonon spectra. We 
have assigned the large peak observed at 33meV to the 
Si (4) mode (for clarity we have numbered the branches 
from highest to lowest in frequency). Although the en- 
ergy of this peak corresponds better to the energies of 
the predicted £4(4) and Si (5) modes, its intensity cor- 
responds much better to that of the Si (4) mode, so we 
are confident that our assignation is correct. 

In Fig.|2|the measured dispersion curves are compared 
with those predicted by the model. The agreement is 
good for the majority of the observed branches. How- 
ever, the model overestimates the energy of branch Si(l) 
by ~ 8 meV (2 THz) at the zone centre and the energy 
of branch Si (4) by ~ lOmeV (2.5 THz) over the whole 
of the Brillouin zone. The relative motions of the atoms 
in these branches are shown in Fig. [3] and it is clear 
that the motion of the copper and oxygen atoms in the 
CuC>2 planes dominates. Both branches change character 
somewhat between the centre and the edge of the Bril- 
louin zone, so the diagrams indicate the motions of the 
atoms at both of these positions in reciprocal space. 

The broad resolution function of the spectrometer 
caused some ambiguity in the assignation of mode sym- 
metries at zone centre. This was due to the possibil- 
ity that we could have accidentally detected out-of-plane 
modes propagating in transverse directions. For instance, 
when we attempted to measure phonons propagating 
along [hhO] it is possible that we could also have picked 
up modes propagating along [00h]. This effect is only 
important at zone centre, so although it could partially 
account for the lowering of the Si(l) mode at zone cen- 
tre, we are confident that it could not have caused us to 
wrongly assign the symmetry of the Si (4) mode. 
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FIG. 3: Anomalous phonon modes. Branches Si(l) and Si (4) 
are shown at the Brillouin zone centre (0 0) and edge (0.5 
0.5 0). The length and thickness of the arrows are approxi- 
mately proportional to the magnitudes of the displacements 
obtained from the common interatomic potential model for 
PrBa2Cu30e+:E (a; = 0). The measured compound [x ~ 0.2) 
has a partial occupancy of the two oxygen sites halfway be- 
tween the Cu sites in the basal plane. 



in PrBa2Cu306.2 the dip at zone centre is not so pro- 
nounced. In contrast to PrBa 2 Cu306.2, the Si (4) branch 
in YBa2Cu306 fits well with the model. Although the 
energies of the other modes are very similar in both com- 
pounds, the Si (4) mode is ~ 12meV (3THz) higher in 
YBa2Cu306 than in PrBa2Cu30g.2 throughout the Bril- 
louin zone. 

The common interatomic potential model assumes that 
the binding mechanism is predominantly ionic, but the 
anomalies we have observed suggest the existence of more 
complex binding mechanisms. It is particularly note- 
worthy that the anomalous phonon branches are charac- 
terised by large displacements of the oxygen atoms whose 
2p orbitals are proposed to hybridise with the Pr 4/ or- 
bitals. Pr-0 hybridisation, in which the Pr-0 bonds are 
partially covalent in character, would result in a modified 
electron distribution around the in-plane O atoms, which 
could change the strength of the Cu-0 bonds and alter 
certain vibrational modes. 



V. CONCLUSION 

We have described inelastic neutron scattering mea- 
surements of the phonon dispersion curves in oxygen- 
deficient single crystal PrBa 2 Cu306+2; (x « 0.2). The re- 
sults have been compared with a model of the lattice dy- 
namics based on a common interatomic potential, which 
is found to overestimate the frequencies of two branches 
by a considerable amount. Strikingly, these branches are 
dominated by motion of the oxygen atoms in the Cu02 
planes: the same atoms that have been proposed to hy- 
bridise with the Pr orbitals. In contrast, the observed 
frequencies of all phonon branches in YBa2Cu30g agree 
well with the model. The frequency shifts observed in 
the two anomalous branches in PrBa2Cu306.2 are inter- 
preted as indirect evidence for hybridisation of the Pr 4/ 
and O 2p orbitals. 



We now compare the dispersion curves we have ob- 
served in PrBa2Cu306. 2 with those previously observed 
in YBa 2 Cu306- 7 In YBa2Cu30g all the branches are well 
described by the common interatomic potential model. 
Therefore the energy shifts of the Si(l) and Si (4) 
branches in PrBa2Cu30g.2 are anomalous. Although the 
Si(l) branch in YBa2Cu306 has a similar shape to that 
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